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This article investigates the generation of non-Gaussianity during inflation. In the context of 
multi-field inflation, we detail a mechanism that can create significant primordial non-Gaussianities 
in the adiabatic mode while preserving the scale invariance of the power spectrum. This mechanism 
is based on the generation of non-Gaussian isocurvature fluctuations which are then transfered 
to the adiabatic modes through a bend in the classical inflaton trajectory. Natural realizations 
involve quartic self-interaction terms for which a full computation can be performed. The expected 
statistical properties of the resulting metric fluctuations are shown to be the superposition of a 
Gaussian and a non-Gaussian contribution of the same variance. The relative weight of these two 
contributions is related to the total bending in field space. We explicit the non-Gaussian probability 
distribution function which appears to be described by a single new parameter. Only two new 
parameters therefore suffice in describing the non-Gaussianity. 
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I. INTRODUCTION 
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The large scale structures of the universe are usually considered to arise from vacuum quantum fluctuations that 
are amplified during a stage of accelerated expansion. In its simplest version, inflation predicts the existence of an 
adiabatic initial fluctuation with Gaussian statistics and an almost scale- invariant spectrum |l| . And it is clear that as 
long as the evolution is linear from the radiation era, non-Gaussianity can arise only from an "initial" non-Gaussianity 
generated during inflation. 

Simple calculations, that we reproduce in the second section of this paper, however show that within single field 



inflationary framework it is not possible to produce primordial non-Gaussian fluctuations if the slow-roll conditions 
Oh, are preserved throughout the inflationary period during which the seeds of the large-scale structures are generated. 
Non-Gaussianity can be generated only if the inflation starts from a non-vacuum initial state |Is| or if there exist 
sharp features in the shape of the potential [|l9j , but it in the latter case it clearly shows up in the density fluctuation 
power spectrum. It has been noticed however that the situation is somewhat changed when more than one light 
scalar field are present during inflation. In this case it is to be noticed first that one generically produces a mixture of 
adiabatic and isocurvature fluctuations |^, [j| that can be uncorrelated or correlated J| | § j . The observational 
consequences of the existence of these two types of modes have started to be considered [l2[ but they clearly 

depend on which type of matter each of the field decays to. Multi-field inflation also opens the door to the generation 
of non-Gaussianity simply because the non-linear couplings can be much stronger in the isocurvature direction than 



in the adiabatic direction 21, |22|, |23|, |24( ]. For instance in models in which a Peccei-Quinn symmetry is broken 
during inflation, it was pointed out by Allen et al. ]l3t that a fourth order derivative term in the effective theory 
leads to non-Gaussianity in the axion density. In the "seed" models including \ 2 [[[fjj axion [|l3| [l5|], Goldstone 
bosons [[l6] and topological defects models, there is a test scalar field that is a pure perturbation and that does 
not contribute to the background energy density; the energy being quadratic in the field, it induces non-Gaussianity 
in the perturbations. In such models however the non-Gaussian features are present in the isocurvature modes only 
and will be observationally relevant only if those modes survive the reheating phase. The phenomenological situation 
is somewhat different however if a transfer of the modes is possible, that is when the fields are coupled |2^, [2(J. 

The aim of this paper is therefore to explore whether in the context of multi field inflation it is possible to generate 
non-Gaussian features while preserving the adiabatic slow-roll type power spectrum and what would be its observa- 
tional signature. We are obviously motivated by the development of Cosmic Microwave Background (CMB) and large 
scale structures observations that offer an opportunity reconstruct the properties of the primordial metric perturba- 
tions. Primordial fluctuations are more directly probed by CMB observations but then the number of modes that 
can be measured is still small. Up to now, no non-Gaussian signature has been detected in either the 4 year COBE 
data |27], ^8) or the 1 year MAXIMA data ^9|. In large-scale structure surveys the number of independent modes 
one can observe is large but the difficulty is that the non-linear gravitational dynamics [|30[ generates non-Gaussian 
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couplings that can shadow the primordial ones |3l[]. One should then rely on a good understanding of the impact of 
the gravitational dynamics on the observations. Cosmic shear surveys might offer one of the most serious opportunity 
to explore such effects in the coming years |52|, [33| . 

We start (Section II) by a general overview of the generation of non-Gaussianity in single and multi-field inflation. 



It will lead us to define a mechanism that can produce such non-Gaussianity. In Section III, we then consider the 
evolution of perturbations in two-field inflation and summarize how isocurvature perturbation can be transferred 
to the adiabatic component when the trajectory in the field space is curved. The isocurvature mode develops non- 
Gaussianity due to self-interaction; we study in Section [f^ the evolution of such a self-interacting field in an expanding 
universe. After having posed the problem for a quantum scalar field in an inflationary background we address this 
issue from a classical point of view. In Section [v| and [vj, we compute the probability distribution function that can 
be obtained in the class of models considered in this article. 

II. OVERVIEW OF THE MECHANISM 

As explained in the introduction, our goal is to design a model that can produce sufficiently large non-Gaussianity 
at least for a band [A, A + AA] of wavelengths observationally relevant, i.e. that corresponds to the large scale 
structure scales. To guide us in this task, we review the properties of non-Gaussianity in single field and multi-field 
inflation models in order to determine whether they fulfill our requirements. 

Let us start by considering a single field, (j>, in slow-roll inflation. The Klein-Gordon equation for its perturbation 
is of the form 

i'cj) + ZH8(f> ~ ^6<p = -V"5d> - 3V'"(S(f>) 2 + ... (1) 
or 

During the slow-roll regime, \H\ <C H 2 and cj> <C 3H(f> so that 

3ff 2 M| ~ V, Hj)~ -V', (2) 
where M4 is the reduced Planck mass, and the slow-roll conditions can be expressed as e <C 1 and \r)\ -C 1 with 

e=M±y, r, = M 2 — . (3) 

In order for the fluctuations of the scalar field to be large enough, one needs the mass of the field to be much smaller 
than H, in which case one gets 5<f> ~ H. To get the correct amplitude for the primordial fluctuations, one needs 
V3/ 2 /(V'Ml) ~ 1CT 5 so that 

H ~ Kr 5 eM 4 . (4) 

The band of wavelengths [A, A -I- AA] corresponding to large scale structures exits the Hubble radius during the e-folds 

N\ = HAt = A In A, (5) 

during which the slow-roll parameter 77 has varied as Ar] ~ (£ — rje)A(f>/M4 : , with £ = Mf (V'" /V). From Eq. (Q), one 
deduces that 

A<jy/M 4 eHAt. (6) 

Now, if the quadratic term in the r.h.s. of Eq (]l|) is dominant over the linear term then V" /V'" < Stfi ~ H and, using 
Eq (Q), one deduces that Ar] > 10 5 r]N\. This will induce a rapid breakdown of the slow-roll inflation. This can be 
understood by the fact that the potential has to be both flat enough for the fluctuations to develop and steep enough 
for the non-linear terms not to be negligible. A way round to this argument is to consider potential with a sharp 



feature 1 19 but in that case the non-Gaussianity is associated with a departure from scale invariance and is located 



in a very small band of wavelengths. 

Thus, one requires at least one auxiliary field. As a second situation let us consider the case in which this auxiliary 
field does not interact with the inflaton so that the potential has the form 

V{<t>, X ) = U{4>) + \m 2 ^ + ... (7) 
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so that the cosmological evolution drives x towards the minimum of its potential x = 0. It follows that the lowest 
order contribution to its energy density perturbation is not given by the standard expression x$X but one has to go 
to quadratic order. Thus, even if the Klein-Gordon equation is linear, the energy density of the field will develop 
non-Gaussianities with a x 2 statistics and with p x ~ H 4 . This energy has to be compared to the contribution of the 
inflaton ~ U ~ H 2 M 2 . It follows that 8p x /8p^, ~ 10~ 5 (m/H) 2 , so that the contributions of the auxiliary field to 
the background dynamics and to the Poisson equation are negligible. \ can develop non-Gaussianities but they are 
not transferred to the inflaton perturbation and gravitational potential. 

The two fields have to interact, a prototypal example being a scattering term of the form p(j> 2 x 2 as proposed in pi[ , 
so that there are a priori two sources of non-Gaussianities: (i) as in the previous case, the density of the field x will 
be quadratic and (ii) due to the coupling, the Klein-Gordon equation for the field 4> will get an source term in the 
r.h.s. of Eq. ([I]). Consider a potential of the form 

1/(0, x) = U(4>) + \^ 2 X 2 + \m\ 2 + ^X 4 , (8) 

so that the effective mass of the auxiliary field is m 2 s = m 2 + p<j> 2 . If initially |0| > |x|, X will roll toward x = where 
it will stay so that <f> drives the inflation and H 2 AI 2 ~ U. For the auxiliary field to develop non negligible fluctuation, 
one needs that m 2 s < H 2 . If p > then this implies that both m 2 and pcf> 2 have to be smaller than H 2 ; since during 
a period of N\ e- folds, using Eq. (||), the variation of the coupling contribution is of order —2p<fiM4eN\, so that it 
can be smaller than H 2 for a number of e-folds 

N x < lQ- 5 H/(p^), (9) 

during which Sx ~ H. When p < 0, the effective mass can be smaller than H due to cancellation between m and 
p(f> 2 but still, the variation of the mass has to be smaller than H so that the condition (||) is a necessary condition 
whatever the sign of p. 

The Klein-Gordon equation ([l]) now has two source terms Si — px4>3x an d S2 = HX 2 &4 1 - The solution will take the 
form 5<f> = 64> u + £0 NG wher e 5<j) H is the solution of the homogeneous equation and <5</> NG a particular solution. As will 
be detailed in Section |VB, one expects that at Hubble scale crossing <5</> NG ~ N\S / H 2 . It follows that for the first 



source term, S<fi NG ~ N\p<f>x 2 /H 2 ~ N\pcf>. The amplitude of the homogeneous (Gaussian) solution is S4> H ~ H so 
that S(j) NG ~ (p()>/H)5(j) H < 10~ 5 5(j) H . The second source term gives rise to a solution 5<j) NG ~ N\px 2 5(j)/ H 2 ~ N\pS(f>; 
the coefficient p cannot be of order unity since otherwise it will give a contribution of order pip 2 to the effective mass 
of x- Thus if p is of order unity either m B s < H so that x ~ H but then it implies that <j)/H < 1, which is not the 
case during inflation or m e s > H but then \ H. It follows that in any situation the non-Gaussian correction is 
negligible compared to the Gaussian contribution. 

Now, let us study the relative magnitude of the different terms entering the Poisson equation. First, U^Scf) ~ eH^M^ 
while the interaction term is of order V^ x xS(f> ~ g 2 <fiH 3 < (10~ 10 /N\)eH 3 M4 and is thus negligible. KxxX 2 
has a contribution of order A% 4 ~ (X10~ 5 )eH 3 M4 which is negligible and a second contribution g 2 <j) 2 x 2 < 
(KT 10 /N x ){(j)/ H)eH 3 M 4 which is also negligible. Note also that the term V x x ~ ™offX 2 < H 4 ~ 10 _5 eiJM 4 , 
so that it is also subdominant. 

It follows that even if x develops non-Gaussianities and is coupled to the inflaton, the requirement that its effective 
mass is smaller than H during a sufficient number of e-folds imply that it is negligible both in the Poisson equation 
and in the inflaton evolution equation. 



As can be seen from the previous example, a quadratic interaction does not fulfill the requirements. We turn to a 
toy model in which the potential, in a neighborhood of the trajectory in field space, takes the form 

V = U^)+m 2 x (0-Mx+^X n - (10) 

n! 

Indeed, this very unusual form with linear coupling tends to show that a particle physics realization of the 
scenario may be difficult to build and the potential ( |Io| ) is just meant to be an effective potential. As in the 
previous example, we have to estimate the order of magnitude of the effect of x both in the Poisson equation 
and in the equation of evolution of the inflaton. Imposing that the drift in the x direction is smaller than in 
the <j> direction during N\ e-folds gives that N\ < (H/ in x ) . The source term in the Klein-Gordon equation 
(|l|) of the inflaton is m 2 x so that it gives rise to a contribution <5</> NG ~ Nxm^/H < 5cf> H which can be of the 
same order as the homogeneous Gaussian solution S<j) H . In the Poisson equation, the only new contribution is 
of order V x< pS(j)x ~ rn^H 2 < (10~ 5 /N\)eH 3 M4 and is thus negligible compared to the contribution arising from 
U. As in the previous example, the self-interaction term will also turn out to be negligible. As a conclusion, 
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FIG. 1: The trajectories of the fields in the plane (4>,x)i once smoothed on a scale R. Before Hubble scale crossing (R < H), 
the trajectories behave quantumly because the quantum fluctuations are active up to scale H and are not smoothed out. After 
the Hubble scale crossing (R > H), the trajectories can be treated as classical trajectories. Note that this transition happens at 
different time for different values of R. The bundle of classical trajectories then evolves in the two dimensional potential and its 
cross section evolves with time. During the bending of the potential valley, the isocurvature modes project on the isocurvature 
modes and there is a transfer because the length of each trajectory is different depending on its position. 

such a potential will give rise to non-Gaussianities sourced by the coupling to the inflaton perturbation. In such 
a scenario, the field x develops its non-Gaussianities due to non-linear evolution and these non-Gaussianities are 
transferred to the inflaton field due to the coupling source term in its Klein-Gordon equation. Note that it does 
not generate directly non-Gaussianities in the gravitational potential since its energy perturbation is always negligible. 

In a neighborhood of the minimum of the potential given by Xo(0) ; the potential can be developed in powers 
of x — Xo(<^) = $X as a n^X n - The first contribution is o^X 2 ! either a 2 > H 2 and the fluctuation of x w ^ be 
negligible or «2 < H 2 and one needs to consider next terms. If the cubic term is dominant then the trajectory Xo{4>) 
is unstable. Now if it is either negligible with the quartic order or of the same order, the latter situation resulting 
in the fine tuning ~ a^H. For the following n > 3 terms, they will be associated with a particular solution of 
magnitude 8<fi NG ~ a n 5x n / H 2 ~ a n H n ~ 4: 5<fi H . The nth order term will have a non negligible contribution to 5<j) only 
if a n ~ iJ 4 ~". Only for n — 4 can a n be a pure number that can be chosen of order unity. This is the only term 
for which the non-Gaussian fluctuations can develop on a long time scale without any tuning to H. This naturality 
arguments lead us to expect the most efficient and only relevant term to consider is an auxiliary field self-interacting 
with a quartic potential, Ax 4 /4!, and coupled to the inflaton p9[ . 

The sketch of the scenario will thus be the following. The auxiliary field develops non-Gaussianities. The bundle 
of field trajectories in the field space {4>, x) is deformed when the potential minimum is bent. According to the sign 
of A, it is focalized (A > 0) or defocalized (A < 0) and the section of the bundle is deformed, as depicted on figure II] 
(this completely analogous to the deformation of the section of bundle of light propagating in a gravitational field)? 
Neighboring trajectories have different lengths so that the fluctuations along each of the trajectories will be slightly 
different. 

To formalize this scenario, we describe in more details the coupling and transfer between adiabatic and isocurvature 
modes. To characterize the non-Gaussianity, one will need to treat the behavior of the quantum fluctuation before 
Hubble scale crossing. But, the trajectories after Hubble scale crossing can be treated classically and the classical 
contribution is expected to dominate over the quantum one if the bending occurs late enough after Hubble scale 
crossing. Note that the case in which A < is less interesting since the trajectory becomes unstable; Xo(^) corresponds 
to a maximum of the potential and we are more in situation like the one of a phase transition. We will however consider 
the case A < while computing the PDF. 



We describe, following mainly ffl, M , in this section the mixing between adiabatic and isocurvature modes in multi- 
field inflation which is one of the major element of our mechanism (see also e.g. for an earlier study). 



III. TWO-FIELD INFLATION: ADIABATIC-ISOCURVATURE TRANSFER 
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The class of multi-field inflationary models can be derived from a Lagrangian for N scalar fields 

R 1 N 

C = ~T^G + 2 S fyW^Vi ~ V (^> • • • . VJv) (11) 

3=1 

where G is the 4-dimensional Newtonian constant, R the Ricci scalar and Greek indices run from to 3. We describe 
the universe by a Friedmann-Lcmaitrc spacetimc with metric 

ds 2 = -dt 2 + a 2 {t)8 ij dx i dx j (12) 

where t is the cosmic time and a the scale factor. We also introduce the conformal time r\ defined by dt = adr). ft 
follows that the Einstein equations take the form 

B° = ^fcti + 2vJ (13) 

N 

H = -47rG£$ (14) 

3=1 

& + 3fT&=-^ (15) 

where a dot refers to a derivation with respect to the cosmic time, t, H = a/a and Vj refers to a derivation with 
respect to <pj. We also set M± 2 = 8itG. 

Let us consider the simplest case in which we have only two scalar fields tpi and ip 2 . We decompose these two fields 
on the direction tangent (cr)and perpendicular (s) to the trajectory as 

^ n \ _ \4,m I -r> 'l SA,m = I sul 

the angle 6 being defined by 



, M(0)( ^ , M{6)=[ ™\ „ (16) 
si \<j>2 J V — sm 8 cos 9 ' 



cos9 = , sin6> = (17) 

\A?l+¥>2 

from which it can be deduced that s is constant along the trajectory. The evolution of the field is obtained by 
combining the two Klein-Gordon equations (pi 



and the Friedmann equation is 



..j , aff^ )+M(fi)(£ ]=0 (18) 



H 2 = ^— [a 2 + 2V(a, s = constant)] . (19) 
3 

The interpretation of the fields a and s as the adiabatic and isocurvature components of 1,(^2) will arise from 
the following study of the perturbations evolution and properties. At linear order, the perturbed metric takes the 
form 0] 

ds 2 = -(1 + 2A)dt 2 + 2a(t)d i Bdx i dt + a 2 (t)[{l - 2C)% + 2d ij E]dx i dx j (20) 
where A, B 7 C and E are four scalar perturbations and we expand the scalar fields as 

<p i (t,x j )=<p i + 6tp i {t,x j ). (21) 
We can introduce two sets of gauge invariant variables. The Newtonian or longitudinal gauge is defined as 

B = E = 0, A = *, C7 = tf, 6tpi=Xi- (22) 
We introduce the perturbation of the inflaton in the flat slicing gauge (C = 0, E = 0) by 

Q l = X* + (23) 
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In Newtonian gauge, the evolution equations in Fourier space reduce to the set 

$ = * 

-4* + + (3H 2 + = --L-(V!Xi + ifiXi + V2X2 + 952X2) 
A 

Xi + ZH± l - -^Xt + VijXj = -2^$ + 4&$. 
By construction, 5s is gauge invariant so that Q s — \s = 5s and from Eq. (p3|), Q a = Xa + cjQ / H . Using that 



t M = -9 t MJ, t M = -9 f MJ-8 2t M : J 



the Klein- Gordon equations (|27|) can be rewritten as 

f 3H 



Xa 

5s 



+MV l] t M 

a 2 



1 

-1 



= 29 J 



(24) 
(25) 

(26) 

(27) 

(28) 



Xa 

Ss 



- 2<$>M 



4$ 



and the Poisson equation (Pq ) takes the form 

A 



47rGa 2 

Defining the mass matrix U of the two fields (a, s) as 



$ = 20<T<5s + axa - <?(Xa - o-$). 



as \ 



= MVi/M, 



and 



= M 



and making use of (pp|) , the system ( |29|) takes the form 

, A 

Xa + 3H±a 



U aa Xa = -2$C/ CT + 4o-$ + 2{6Ss)- - 2-U a 5s 



Ss + 3HSs 



U ss )Ss 



6 A 



a 2irGa 2 



The comoving curvature perturbation is given by 



U H H 

H = C A —(0i^i + 1P2S1P2) = C + —Scr 

p + f 



a a 



from which it is deduced that 



TT A TJ 
U= -^^$ + 2— 05s. 
Ha 2 a 



(29) 

(30) 

(31) 

(32) 
(33) 

(34) 
(35) 



It follows from this analysis that as long as 6 = 0, 5s evolves as a test scalar field evolving in an unperturbed 
Friedmann-Lemaitre universe and does not affect the evolution of the gravitational perturbations; $ couples solely 
to the fluctuation of a in Eq. (|30|). It is recovered that the entropy remains constant on super-Hubble scales as 
long as there is no mixing. 5s will transfer energy to the gravitational potential only when 9 7^ and the comoving 
curvature can be affected significantly even on super-Hubble scales if the fields follow a curved trajectory. Note that 
this mechanism is effective only during the inflationary period. 

The "surviving" isocurvature perturbations (all of them in the case where there is no bending) needs to be taken 
into account. At the end of the inflationary period, both the inflaton and auxiliary field will decay into particles and 
radiation. To set the initial conditions in the radiation era, one needs to specify in details these decays. When there 
is no bending of the trajectory, this is the only mechanism through which an imprint of the isocurvature modes can 
survive. According to the scenario, the initial perturbations in the cosmic fluids at the beginning of the radiation era 
are a mixture of adiabatic and isocurvature modes that can be correlated (see e.g. j(| |^]). In full generality, we will 
have to consider both contributions but we focus in the following on the "pure" case where the fluctuations at the 
end of inflation are strictly adiabatic, which can occur if the two fields decay identically. 
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IV. GENERATION OF NON-GAUSSIANITY 



The goal of this section is to estimate the magnitude of the non-Gaussianity developed by the isocurvature mode 
during a phase of de Sitter inflation. The frame in which such calculations should be performed is the quantum field 
theory for a coupled scalar field. As we are interested in the emergence of weak non-Gaussian features, through for 
instance the emergence of non-zero connected part of high order correlation function, a perturbation theory approach 
should be applicable in principle. 

We point out however that a series of technical or conceptual problems emerge in this physical situation for which 
there seems to exist no known solution. We first present in which way the computation we would like to do is 
affected by those problems but, as their resolution goes far beyond the estimate we would like to obtain, we finally 
turn to a classical treatment of the field behavior. This provides a solid enough ground to estimate the amounts 
of non-Gaussianities at scales that remain super-Hubble a time long enough before the adiabatic-isocurvature mode 
transfer. 



A. The quantum level calculation 

As seen from the previous investigation, the entropy field Ss is decoupled from the gravitational perturbations and 
can be considered as a test field evolving in an homogeneous and isotropic cosmological universe. To characterize 
the statistical properties of such a self-interacting field, one would ideally like to compute its high-order correlation 
functions. The inflationary phase can be described by a de Sitter spacetime |34| in flat spatial section slicing 

ds 2 = (-drf + S ij dx i dx j ) (36) 

which is conformal to half of the Minkowski spacetime. The conformal time is related to the cosmological time by 

V = ~e- Ht (37) 

and runs from — oo to 0, the limit 77 — s- — representing the "infinite future". The de Sitter spacetime can be viewed 
as a four dimensional hyperboloid embedded in a five dimensional Minkowski spacetime 

- M 2 + (zi) 2 + {X2? + (x 3 f + (.x 4 ) 2 = l/H 2 . (38) 

The invariance of this surface under five dimensional Lorentz transformations implies that the de Sitter space enjoys 
a 10 parameter group of isometries known as the de Sitter group 0(4, 1). It is usual to define the de Sitter length 
function z(x, y) by H 2 (x — j/) 2 /2 = 1 — z(x, y) where x and y are the five dimensional coordinates of two points on the 
hyperboloid (pq). It follows that z(x,y) — H 2 xy and the two points are timelike or spacelike separated respectively 
when z > 1 and z < 1 

For a minimally coupled free quantum field of mass m, due to the spatial translation invariance, the solution can 
be decomposed in plane waves as 

/d 3 k r ^ -m n 

[vo(k )V )b k e*-* + vZ(k, V $le- ik '*] (39) 

where we have introduced v = a 5s, a hat referring to an operator. In this Heisenberg picture, the field has become 
a time-dependent operator expanded in terms of time-independent creation and annihilation operators satisfying the 
usual commutation relations [&k, = (27r) 3 5^ 3 ^(k — k'). We can then define the free vacuum state by the requirement 

6 k |0) = for all k. (40) 

As it is standard while working in curved space |]35| , the definition of the vacuum state suffers from some arbitrariness 
since it depends on the choice of the set of modes Vo(k, if). They satisfy the evolution equation 

t *_J._^Wo, (41) 



77^ rf- 



the general solution of which is given by ^J-K-q/A c\Hy (kij) + C2HI (JtT)) with | C2 1 2 — |ci | 2 = 1, where Hi ' and 
are the Hankel functions of first and second kind and with v 2 = 9/4 — m 2 / H 2 . Among this family of solutions, 
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it is natural to choose the one enjoying the de Sitter symmetry and the same short distance behavior than in flat 
spacetime. This leads to 

v Q {k,r 1 ) = ^^Hi 2 \kri). (42) 

This uniquely defines a de Sitter invariant vacuum state referred to as the Bunch-Davies state vacuum |3^]. In the 
massless limit, the solution ([Il|) reduces to 

i \ c~ ikl > 



•"'"I'-sjJw (43) 

Having determined the free field solutions, one can then aim to express the iV-point correlation functions of the 
interacting field, Ssi, in terms of those of the free scalar field. For instance, on the example of a Ss 4 theory, the 
connected part of the 4-point correlator of the interacting field will reduce, at lowest order, to 

(0|T& l (x 1 ,r ?1 )...& i (x4,r ?4 )|0) =i- / (0| T&(xi, r?i)&(x 2 , 77 2 )&(x 3 , r? 3 )&(x 4 , m)& (x, rj) |0) /=^d 4 xd7?. (44) 

Using the Wick theorem and keeping only the connected part leads to 

(0|T& i (xi,77i)...& i (x 4 ,77 4 )|0) =i- / ]^[G(x,77;x,,77 4 )y^d 3 xd?7. (45) 

' "* i=i 

G(x, 77; x', 7/) is the free propagator, that is the time ordered product of two free fields in the free vacuum 

iG(x,r,;x',V) = a( } * ( ,) / 7S5 M*. vK(k,v')e^- Ax 0(-A v ) + v (k, 7 7 ')uS(fc,7 7 )e ik - Ax 0(A7 7 )] (46) 

with Ax = x' — x and A?7 = 77' — 77 and 9 being the Heavyside function. After integration over angles the free 
propagator can be computed [B7| to be 



1 G(x,r ? ;x,r ? ) = i ^^ Ax2 _ A?72 +l + y o — cos(fcAa;)e~ 



ifc|Ar) 



(47) 



While trying to compute the 4-point correlator ( f46| ) with the latter expression of the free field propagator ( Jl7j ) , one 
has to face the existence of two infrared (IR) divergences. Note that at coinciding points, there is a UV divergence 
that can be regularized by taking into account that inflation started at a given initial time |3(| . 

The IR logarithmic divergence at k — was first exhibited by Ford and Parker (3^]. Allen and Folacci |3j| showed 
that this divergence arises from the incorrect assumption of de Sitter invariance for the vacuum and from the existence 
of a zero mode, i.e. the action is invariant under transformations of the form <f> — > <f> + constant. It is well known that 
an expansion in terms of creation and annihilation operators as in Eq. ( |39|) is inadequate for the zero modes |40| in 
the same way as the expansion in terms of creation and annihilation operators for the standard quantum harmonic 



oscillator breaks down when the frequency is going to zero. To be slightly more precise, let us recall 1 41 the properties 
of a free massless scalar field leaving on a 3-torus of volume V. Beside the standard plane wave solutions associated to 
the creation and annihilation operators, 6k and 6^ obtained for k and describing harmonic oscillator of frequency 
|k|, a complete set of solutions requires to consider the solution (xq +pot)/vV obtained for k = and describing the 
classical solution for a free particle. The position and momentum operators Xq and po satisfy the commutation relation 
[icoiPo] — i- A vacuum state can then be defined by imposing p |0) = and |0) = 0. This state is the product of 
Fock and a Hilbert space corresponding respectively to the oscillators and the free particle (see e.g. f|l, 42, 43|, fS)) 



and will not be normalizable. For a flat space in more than two dimensions, the continuum limit exists because the 
contribution of the zero mode is of zero measure, the volume of the phase space cancelling the divergence but the 
effect remains in de Sitter space whatever its dimension. 

This divergence led some authors fl3L [f5| to define other vacua with less symmetry than the de Sitter group but 
with a well defined propagator. For instance, in the closed spatial section slicing a natural choice is the 0(4) invariant 
vacuum that is symmetric under rotations of the constant time hypersurfaces. The modes are discrete and the 
IR divergence is avoided by choosing a set of modes with a different solution for k = 0. Kirsten and Garriga Q| 
proposed the construction of an acceptable de Sitter invariant vacuum in which the zero mode is well treated. The 
case of de Sitter space with static spatial sections slicing was considered by Polarski E5I . In the case of flat spatial 
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sections slicing, which we are most interested in, the IR problem can be regularized by working on a torus | p7| with 
—H~ 1 /2 < x 1 < H~ 1 /2 which is equivalent to set an infrared cut-off. With such a regularization, this yields the 
result 



iG(x, rj; x', rf) = |- ^_!^__ _ I m H ^ Ax 2 - A V 2 ) + constant ) . (48) 



A second IR divergence arises at late time, i.e. when 77 — > since the integral (|4(]) with the regularized propagator 
(|48f) diverges due to the volume factor l/iJ 4 r/ 4 . The former resolution of the IR divergence in k = cannot resolve 



this late time IR divergence. Tsamis and Woodard |37 pointed out that such correlation functions suffers from a 
series of flaws: (1) they are not finite even at lowest order, this problem becoming worth as the number of vertices 
grows, (2) it is not purely imaginary so that it implies a tree-order breakdown of unitarity. The physical origin seems 
to be the redshifting that drives all physical momenta toward zero when 77 — > 0, making the overlap between plane 
waves very strong. It is to be noted that the use of the regularized propagator given in Eq. ( f48| ) does not cure the 
problem nor the introduction of a late time cut-off (the existence of which could be associated with the reheating 
time). 

The resolution of these fundamental problems goes far beyond the estimates we want to obtain on the effects of 
nonlinear couplings. We thus adopt a simpler approach assuming that at scales that exceed the Hubble size the field 
value trajectories are classical (e.g. deterministic) and encoded in the potential shape. Because the trajectories might 
have a non-trivial dependence with the initial field values set up at Hubble scale crossing, non-Gaussianities can be 
induced during that period. If the time between Hubble crossing and the bending of the trajectory is long enough the 
non-Gaussianities of such classical origin are going to exceeds those a priori present in the initial value distribution 
as tit emerges from the quantum process. 



B. The classical limit 



To implement this idea, we consider the field & at a large enough scale. It amounts to applying to the evolution 
equation of the field a filtering procedure at a fixed (comoving) scale, R. In the following we note Ss the filtered field 
(what exactly is the function it has been convolved with is not important). At the time the Hubble size has shrunk 
below the smoothing length, the trajectory of the field becomes classical. Its evolution equation is simply, 

Ss + 3HSs = S{5s). (49) 

It derives from the real space Klein-Gordon equation applied for Ss [e.g. Eq. (p3[)] where the Laplacian term has 
been dropped on super-Hubble scales. Assuming that the trajectory of the mean field value in a super-Hubble patch 
is insensitive to the small scale fluctuations we can replace S_(Ss) by S(&)_ In the absence of source terms 
in Eq. (f49|), Ss is simply constant (e.g. trajectories are parallel lines on Fig. |l]) and its value is given by its initial 
value set up at Hubble crossing. The initial value is given by a sum of Gaussian distributed values and non-Gaussian 
corrections induced by the non-linear couplings during the sub-Hubble evolution of the field, 

& init . =& (0) +& (ng) . (50) 

The resolution of Eq. ([l9j) requires in general the knowledge of the source term. It can be solved however pertur- 
batively if one assumes that Ss can be expanded in terms of the coupling constant entering the source term, 

&(*) = * (0) + •■• (5!) 

The evolution equation for is trivial to obtain. It reads, 

Ss (1) +3H5s (1) =S(5s {0) ) (52) 



the solution of which is 



which can be rewritten as 



Ss^(t)=&^ + (t-t K )^p-., 



S(5s {0) ) 



(53) 



&W(t) =&<»«> +JV*^_i, (54) 
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where Nr is the number of e-folds from the Hubble size crossing, tjj, to the time t. If the latter is large enough one 
expects the second term of this equation to dominate the first, e.g. that non-Gaussian effects induced during the 
classical evolution dominates those built during the stage of the quantum evolution. We will make this hypothesis in 
the following. 

It clearly implies, as we anticipated in part II, that the amplitude of the first order corrective term is given by the 
amplitude of the source term S(5s^) divided by H 2 times the number of e-fold during which the non linear couplings 
are active. 

The existence of such corrective terms are obviously of importance for the statistical properties of the field. If the 
source term contains nonlinear couplings the field 5s is no more Gaussian. Such properties should be exhibited in the 
high order correlation function of the field. Exploring these consequences is the aim of the next section. 



V. ISOCURVATURE MODES PROBABILITY DISTRIBUTION FUNCTION 



As seen from the previous analysis, nothing prevents the isocurvature modes to develop non-Gaussian properties 
on super-Hubble scales. As stressed before, if the time between Hubble radius crossing and the exchange of modes is 
large enough, the properties of those modes will be determined by their stochastic evolution at super-Hubble scale, 
not so much by the quantum state with which the field modes reach super-Hubble scales. In other words we expect 
the high order correlation functions present in the quantum field to be finally superseeded by the ones induced by the 
subsequent stochastic couplings. 

The aim of this section is then to characterize the way the non-Gaussian features in the isocurvature modes are 
built up. There are obviously many ways of characterizing non-Gaussian features in a stochastic field. The simplest 
approach is to consider the shape of the one-point probability distribution of the local field value. This series depends 
obviously in the type of couplings one has. For the reasons previously detailed, we will consider only the case of 
quartic couplings. 

In the first part of this section (§ VA) we show how it is possible to compute any of such cumulants at leading 
order in the coupling constant (that would correspond to tree-order calculation in a quantum field formulation). It 
finally leads to the expression of the generating function of the one-point cumulants (§ |VB|). 



To get insights into the physical interpretation of the previous results we then present (§ VC) the derivation of a 
quantity more directly related to observations: the one-point PDF of the local field value. Not surprisingly we will 
see that the rare event tails of those distribution differs from those expected for a Gaussian distribution. Depending 
on the sign of the quartic coupling, one expects an excess or a deficit of rare values. 



A. The expression of the fourth order cumulant 

Let us focus on the case of a stochastic field 5s self-interacting with a potential, 



V(&) = ±5s i (55) 

and we assume that the coupling constant A is small compared to unity. In the 5s expansion, Eq. (|l]), it means that 
is of the order of A™. We then consider the evolution equation ( ^9|) with 

S(Ss)=-^5s 3 . (56) 

At this stage it should be noted that when using the 5s field in the source term, one neglects the effects of the small 
scale fluctuations in the field trajectory. This is of no consequence for tree order calculations, however, no loop terms, 
that involve arbitrarily small scales fluctuations, can be reliably computed from this approach. If one wants to do that 
one should solve the full quantum problem. In this study we allow ourselves to use this simplification but forbid 
ourselves to compute correlation properties beyond tree order. 

As said before in the absence of coupling the free field solution Ss^ is time independent. This is not the case of 
higher order terms. For instance from Eq. (154), we know that for this coupling term, 



& (i) --4^ 



18 H 2 



(57) 



a term which induces non-Gaussian properties in the field. In case of a quartic coupling the value of 5s remains 
symmetric distributed so that the third order cumulant is always zero. The first non-vanishing high order cumulant 
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FIG. 2: Diagrammatic representation of the fourth order cumulant in a perturbation theory approach. Lines correspond to 
connected pair points in ensemble averages of product of Gaussian variables as an application of the Wick theorem. End points 
are taken at zero order in the coupling constant (they are linear in the initial Gaussian field). Points at first order in the 
coupling constant are cubic in the field, at second order they are quintic, etc. The first diagram corresponds to the tree order 
term. It involves only one vertex. The other two diagrams correspond to loop corrections. 




FIG. 3: Diagrammatic representation of the sixth order cumulant at leading order. Two types of diagrams appear, one with 
two three-leg vertices, one with one five-leg vertex. 



is then the fourth one. Its value can be computed at leading order (in the coupling constant) from the expression of 
Indeed the fourth order cumulant is formally given by, 



<& 4 >cM<F)-3<&y 



(58) 



replacing Ss by its expansion in the coupling constant, Eq. (|51|). The computation of these moments up to linear 
order in A gives, 



As 



(1) 



3(Ss [1 >Ss 



As 



(0) 



(59) 



This first two terms of this expression cancel because Ss^ ' obeys a Gaussian statistics. One can finally check that the 



latter expression gives the connected part of ( Ss^ 
given by, 



so that, at leading order, the fourth order cumulant is 



(& 4 ) c = 4(&« 



which can be easily computed from the expression of , so that, 



<& 4 > c 



4A N R 



(60) 



(61) 



It is clear that at leading order the expression of the fourth order cumulant involves only the expression of Ser"' . 
All higher order terms in the expression of Ss contribute only at higher order in A to this cumulant, and they would 
correspond to loop term corrections, e.g. Fig. ^[ 

However if one wishes to compute the expression of higher order cumulants, such as the sixth order one, the higher 
order terms in the expansion of Ss will then contribute (still at tree order!), see Fig. In the next subsection we 
explore in more details the mathematical structure of the cumulant hierarchy. 



B. The structure of the field cumulants 



To get insights in how the leading order expression can be computed, it is useful to consider one where terms of 
various order in the expansion of & are mixed together. 



12 



For instance the sixth order cumulant is, at leading order in A, given by 



<* 6 > c 



30 



(i) 



<fe 



(0) 



6 6s {2) 



(0) 



(62) 



since Ss^ scales like A 2 &^ . The factors that appear in these expressions correspond to the number of each of 

such terms that appear in those expansions. 

In general the tree order term of any cumulant is determined by the consequences of the Wick theorem applied to 
Ss^\ It will therefore be of the form, 



?™>c= yi (n> (ps) •■•* (p " ) 

decompositions \i— 1 / 



where all the decompositions are such that 



]T(2 Pl + 1) = 2(n - 1) 



(63) 



(64) 



so that it is possible to connect all terms together and no loop can be built [2(n — 1) is two times the number of lines 
required to connect n points] . 

In order to manipulate dimensionless numbers we introduce the (time dependent) vertices, 



t>+i = (2p + 1) 



6s (p) 



2p+l ■ 



For instance 



At 



Nr 
3H 2 ' 



Then the high order cumulants can be rewritten as 



<a= E n 

decompositions Li=l 



2pi+l 



n-1 



(65) 



(66) 



(67) 



that appear to be identical to tree sums in which the p-leg vertices are v p . The general computation of the v p is still 
a difficult task since it requires the resolution of the evolution equation for 6s. Actually if the time after Hubble scale 
crossing is large enough then the first term of the equation is finally negligible so that the evolution equation actually 
reads, 



3HSs(t) = -^ls(t) 3 . 
Unlike the full evolution equation this equation has a simple solution given by, 



[ss 2 (t)V = 



A (t — t H ) 



36 H 



the time dependence of which can be rewritten as a function 1/3 (t) , 



&(*) 



1 - Mt) 



1 2 



(68) 



(69) 



(70) 



/3 



where &^ ^ is the (time independent) value of 5s at Hubble scale crossing. It is to be noted that the direct use of 
this relation to compute the high order moments of Ss is unjustified because it automatically introduces loop terms 
in the computations. As these terms cannot be reliably computed we restrict ourselves in the following to tree order 
terms. We can also note that if ^3 is positive there is a singularity at finite distance which would make all moment 



13 



values infinite in this case. This pathological behavior is due to the fact that on rare occasions, the field value Ss 
moves at arbitrarily large distance from the origin in a potential which is unbounded from below. This instability is 
due to the description of the potential we use around the trajectory and is not necessarily physical. Once again if 
one limits ourself to tree order computations this issue is automatically solved; all cumulants are finite are tree order. 
This calculation rules provides us with a natural regularization scheme. 

A convenient way to describe what we have obtained is to write the fourth order cumulant, in units of the second 
order moment, 

±-k = 4va. (71) 
The sixth order cumulant can similarly be expressed as a function of ^3 and v§ 



2\ t > 



(30i/5+6i/&). (72) 



In general the cumulant of order 2n scales like (<fe 2 ) to the power In — 1 with a ratio given by a sum of product of 
vertices. 

It is useful to define the vertex generating function, 



Q(t,r) = J2M^ (73) 
v 

which can be straightforwardly related to the expression of 5s in Eq. ( |70| ) when Ss^ is replaced by r so that 

G(t, t) = T (74) 

This expression provides the values of all the vertices that can be expressed in terms of 1/3 (i), at least for large values 
of Nr. We emphasize that this technique can be applied to any potential shape, although it is probably not always 
possible to find a close form for the generating function in all cases. 



C. The PDF calculation 

1. General formalism 

The computation of the one-point PDF of the isocurvature fluctuations relies on the use of its cumulant generating 
function, (p(y) defined below. All cumulants can obviously be obtained, order by order, from the vertex generating 
function, Q, we have just obtained. But it is actually possible to take advantage of their tree structure to compute 
the whole cumulant generating function at once. This latter is defined as, 



(Ss 2n )c {-yf 

This latter, because of the tree structure we are dealing with, is obtained from a Legendre transform of the cumulant 



V(y) = l^ , 2 2n-l 2nl (75) 



generating function 1 30 



where r(y) is solution of 



<p{y) = yG(T{y)) + \r 2 (y) (76) 



r(y) = -y^. (77) 



The derivation of the equation system is too long to be recalled here (see Q f° r details). It is however worth noting 
that 



(78) 
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FIG. 4: Integration path for Eq. ( |79| ) in the y-complex plane (thin solid lines) for v% H 2 = 0.3 and 5s/a s = 0.3 and 1. The two 
half bold lines on the real axis represent the location singularities for ip{y). For large values of <j> the integration path is pushed 
along the singular part (dashed line). 




FIG. 5: Shape of the one-point PDF of Ss for z/3 H 2 = 0.3 (dot-dashed line) or vz H 2 — —0.3 (solid line) compared to a Gaussian 
distribution (dashed line). 



The one-point probability distribution function is then given by the inverse Laplace transform of ip{y) 

' 3 dy 



P(&) 



O • 2 eX P 



y{y) Ss y 

rr 2 ' „ 2 



(79) 



where a s is the variance of Ss. The global properties of P(Ss) can be derived from the properties of the cumulant 
generating function. 

From the expression (174) one can obtain the expression of r in @, 



r(l - is 3 t 2 /3) 



3/2 



-y 



(80) 



from which the expression of tp(y) can be explicitly computed. The properties of tp(y) in the complex plane are going 
to depend on those of r. 

The shape of the PDF can be calculated from a numerical integration in the complex plane. In practice to complete 
such a numerical integration one must choose an adequate path in the y plane to make the integral convergent. This 
is achieved in imposing that the quantity in the exponent remains real along the trajectory which can be obtained if 
the trajectory crosses the real axis at the saddle point position, y 8 , defined by 



dy 



y=y 3 HP 1 



(81) 



which with Eq. ( ]?8\) is given by the value of t s for which Q{t s ) = Ss and then y s = —t s /Q'{t s ). Example of such 
integration path are presented on Fig. [|. 
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The resulting PDFs are shown on Fig. [| They clearly exhibit non-Gaussian features (the Gaussian case is shown 
as a dashed line). The magnitude of these features depend on the value and sign of V3, that is indeed related to the 
value and sign of the coupling constant A. In the following we explore in some details the behavior of the PDF in 
different cases. 



VI. PROPERTIES OF THE ISOCURVATURE MODE PDF 



The behavior of the PDF for small values of Ss 



For small values of Ss/a s it is possible to derive an explicit expression for the one-point PDF based on a saddle 
point approximation. In Eq. ([79]) one can expand (fi(y) around the saddle point defined previously, Eq. (|8~l|). 
The resulting formal expression of the PDF then reads, 



P{Ss)dSs 



dSs 



which can be calculated from the explicit expression of tp, 



exp 



P(S,)dSs=J- 



1 — Ss 1/3 



(3+<r^ 3 ) s 



exp 



2<7„ 



3& 2 



(6 + 2& 4 vt)<r, 



d^ 
as ' 



(82) 



(83) 



This result is valid as long as Ss 2 U3 is small compared to unity. For larger values of Ss the behavior of the PDF 
depends crucially on the sign the A. 

This formula provides actually a very good description of the overall PDF shape when A is positive (since excursion 
to rare events are anyway not permitted). 



B. The rare event tails in case of negative A 



When 1/3 is positive, that is when A is negative, the shape of the potential is such that it favors the rare event 
tails. The behavior of the probability distribution function in the rare event tails depends in particular on the analytic 
properties of ip(y) in the complex plane. For positive values of 1/3 it can be easily checked from Eq. (|8^) that r(y), and 
consequently (p(y) is non-analytic on the real axis with two symmetric singularities at finite distances of the origin, 

t c = 4a/^, (84) 
2 V V3 

(e = ±1) corresponding to values of r where its derivative with respect to y is diverging. In the vicinity of this point 
it is easy to expand first the expression of y as a function of r, 

y = ey c + e (t - t c ) 2 (85) 

with 

y c = --^— (86) 



and then the expression of f(y) reads [taking advantage of Eq. [78| 

ip(y) = ip c + r c (ey - y c ) - a c (e y - y c ) 3/2 + ... (87) 



with 



= -^r- (88) 

16^3 

' (89) 



(90) 



V3 
16 



9\/3> 



1/4 
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Those singularities induce exponential tails in the PDF of & (a very complete mathematical investigation of such 
cases can be found in jf?}). The shape of the tails are of the form exp(— \5s\/6s c /cr s 2 ) where Ss c is given by the inverse 
value of y at the singularity, & c = 16 A /7/3~/9. 

To be more precise the expression (|83| ) ceases to be valid for |&| > Ss m where Ss m = Q{t c )- For larger values of 
& the saddle point position in y is pushed towards one of the singularities, ±j/ c . The behavior of the PDF will then 
be dominated by the behavior of <p(y) around this point, Eq. (|87|), in which a regular part (p c + r c (ey — y c ) and a 
singular part a c (ey — y c ) 3 ^ 2 appear. One can then make the integration path in y in Eq. ( f79| ) running along the real 
axis (both ways) crossing it at y = y c so that it can be described by the real variable u varying from to oo with, 



V = e y c + e u e 



(91) 



where the sign is changing according to whether y is above or under the real axis. Expanding the singular part in the 
exponential one gets, 



P(&) 



2tti (7, 



. u 3/2 e ±3W2 exp 



-ip c + ey c Ss - (e&s - r c )u 



(92) 



The singular value in y to consider (that is the value of e in the previous equation) depends on the sign of 5s (the 
positive tail corresponds to negative value of y) so that finally one gets, 



P(5s) = 



(\Ss\ - r c ) 5/2 exp ^- 



ip s + \y s 5s\ 



r(-3/2) 

which, for the parameters describing the singular behavior of </?, gives, 



P(Ss) 



4 a, 



3v37r v. 



1/4 



Ss\ 



1 



-5/2 



exp 



16^3 a s 



9& 



16 v /T^cr s 2 



(93) 



(94) 



It is clear for this expression that the rare event tails are very different from a Gaussian distributed variable. 



C. Bounding values in case of a positive A 



When is negative, that is when A is positive, we are in the opposite case. We expect that the rare event tails 
to be chopped out. Actually in this case it is easy to see that there is no singularity on the real axis. It implies that 
the integration contour can be moved from the imaginary axis to the left or the right to an arbitrarily large distance. 
When |y| — > oo, we have r ~ |y| 1 ^ 4 (— S/v^) 3 / 8 which implies that, 



<p(y) 



3 

-"3 



1/2 



\y\ 



(95) 



for large values of y. As a consequence for large enough values of <f> the integral simply vanishes away. The values of 
<j) are therefore are bounded by ±-\/— 3/1/3. The shape of the PDF near the bounding values can also be computed 
explicitly from the behavior of tp for large values of y, 



v(y) = - — 

Then the expression of the integral ((79|) reads, 

dy 



1/2 



33/4 



9 

i L>3 



P{Ss) = 



27ri a J 2 



exp 



with 



(V^Wi-mM-aW 2 - 9/(8*3) 



33/4 



(-* 3 ) 3 / 4 ' 

A simple change of variable, i 1 / 2 = ay 1 ! 2 j o 2 , shows that it can be written, 

2 / n \ r+i°° 



P(_ &) = fVexp(-^) / |Up(- t V* 
a 2 \ 8is 3 a s 2 J J_ ioo 2m \ 



+ zt 



(96) 



(97) 



(98) 



(99) 
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with 



^ 2 (V /Z 3M-|&I) 



(100) 



The behavior of the PDF near the bounding values are determined by the small values of z. The expression of the 
PDF can be obtained in this regime by a saddle point approximation, similar to Eq. (B3), which leads to, 



P{Ss) 



1 



a 2 20T 



exp 



1 



4 z 8v 3 <t s 



(101) 



which reexpressed in terms of 6s gives the behavior of the PDF near the bounds, 



P(&) 



2a/7TCT s 



exp 



MM) 3 / 2 



4<7 s 2 (v^37^-|&|) 



(102) 



The PDF is found to go continuously to zero at 5s = ±y— 3/^3 positions. 



VII. CONCLUSIONS 



In this paper we have explored the possibility of generating significant non-Gaussian initial metric perturbations in 
the context of inflationary cosmology while preserving a power spectrum of slow roll type adiabatic fluctuations. We 
found that the only viable mechanism is through a multiple field inflation where transverse (e.g. isocurvature) modes 
developed non Gaussian properties that can be subsequently transferred to the adiabatic fluctuations if the classical 
field trajectory is bent. 

We have pointed out that quartic type coupling in the transverse modes is the most natural type of couplings for 
providing non-Gaussianities in the sense that in this case no fine-tuning in the value of the coupling constant has to 
be invoked. We stress that in the context of such a mechanism, unlike any others, the amount of non-Gaussianities 
that can be fuelled in the adiabatic fluctuations can be almost arbitrarily large. 

We have examined in more details the case where the isocurvature mode generation (that is when they reach the 
Hubble size during the inflationary period) and the adiabatic-isocurvature mode mixing happen at very different time. 
The reason we consider this case is two-fold. First it is somewhat pedagogical since it shows that these two stages do 
not have to be concomitant. Second it implies that the non-Gaussian properties of the isocurvature modes developed 
mainly during the time they live at super-Hubble scales. It makes their computation much more simple since we can 
avoid a full treatment of the nonlinear field evolution at a quantum level (and it turns out that such a computation is 
not straightforward at all!). For modes living at super-Hubble scales we allowed ourselves to view the field evolution 
as the one of a classical stochastic field. In this case it is then possible to pursue the calculations to completion in the 
sense that it is possible to derive their whole set of correlation properties. 

In particular we have been able to derive the one-point field cumulants at tree order in the weak coupling limit in 
a consistent way and finally build up the one-point probability distribution function of the field value. In such class 
of models, the statistical properties of the curvature perturbation are described by the superposition of a Gaussian 
and a non-Gaussian contributions with a relative weight proportional to the bending angle, A6, of the trajectory in 
field space during slow- roll. The non-Gaussian component is fully characterized by a single parameter, 1/3, related to 
the reduced fourth order connected cumulant and has the same variance as the Gaussian contribution. For practical 
purposes, we emphasize that its PDF is well approximated by Eq. (|83|) that reproduces accurately the numerically 
computed PDF within our approximation scheme. Thus, all the statistical properties can be encapsuled in two 
parameters. These results give some insights on what type of non-Gaussian features can appear in future large-scale 
structure or CMB surveys while assuming the inflationary prejudice. Note that the kind of non-Gaussianity described 
here departs from that generated by the non-linear gravitational dynamics, in particular it has no skewness. There 
is still however some ways between these results and their observational consequences. How non- Gaussian properties 
that are present at super-Hubble scales are transferred for instance to the CMB anisotropies at sub-Hubble scales is 
not totally straightforward. Whether such effects could be actually observed when observational aspects are taken 
into account demands in-depth analysis. 
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